Phase transitions in the early and the present 
Universe: from the big bang to heavy ion collisions^ 

D. Boyanovsky 

Department of Physics and Astronomy, University of Pittsburgh, 
Pittsburgh PA. 15260, U.S.A 

February 1, 2008 



Abstract 

In these lectures I discuss cosmological phase transitions with the goal of establishing the possibility 
of observational consequences. I argue that the only phase transition amenable of experimental study 
within the foreseeable future is that predicted by QCD and discuss some of the potential observational 
cosmological consequences associated with this phase transition(s). I describe the experimental effort 
to study the QCD phase transition(s) at RHIC and SPS and summarize some of the recent experi- 
mental results. The possibility of novel phases of QCD in the core of pulsars is discussed along with 
the suggested observational consequences. A brief review of standard big bang cosmology as well as 
the astrophysics of compact stars sets the stage for understanding the observational cosmological and 
astrophysical consequences of phase transitions in the standard model. 
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1 Prologue: Cosmological Phase Transitions, Theory vs. Obser- 
vations 

The theme of this School is Cosmological Phase Transitions, Theory and Observation. 

A wealth of cosmological data is providing confirmation of sound theoretical ideas in early Universe 
cosmology. Measurements of temperature anisotropy in the CMB by satellite, balloon borne and earth 
based observations as well as measurements of the acceleration of the expansion with supernovae Type la 
searches and precision measurements of light element abundances (for recent reviews see Q) provide an 
impressive body of complementary high quality data. 

While the unprecedented quantity and quality of cosmological data seems to validate the main ideas of 
early Universe cosmology, the observational consequences of cosmological phase transitions are still rather 
indirect. Although there is more theory than direct observation of aspects of cosmological phase transitions, 
I will argue in these lectures that current and future accelerator experiments along with observations of the 
properties of pulsars will be opening a window to the early Universe at a time scale of about 10~^ seconds 
after the Big Bang. 

To set the stage for the description of the experimental effort to probe the last phase transition predicted 
by the standard model of particle physics, and to understand its potential cosmological impact, 1 review the 
standard hot big bang model and the astrophysics of compact stars, emphasizing the different scales and the 
observations associated with relevant phenomena. 



2 The Standard Hot Big Bang 
2.1 Ingredients 

The standard "Hot Big Bang" theory of Early Universe cosmology is based on a wealth of observations and 
rests upon the following pillars: 

• Homogeneity and isotropy: on large scales > lOOMpc the Universe looks homogeneous and isotropic. 
This is confirmed by galactic surveys of large scale structure and by the homogeneity and isotropy of 
the Cosmic Microwave Background (CMB). 

• The Hubble expansion: objects that are separated by a (comoving) distance d recede from each other 
with a velocity v = H d with H the Hubble parameter (or Hubble constant) whose value today is 
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Hq ~ 65Km/s/Mpc. The Hubble law of expansion determines the size of the causal horizon, objects 
separated by a comoving distance 



h = 



SOOO/i-^Mpc 
H 



(2.1) 



(2.2) 



lOOKm/s/Mpc 



recede from each other at the speed of light and are therefore causally disconnected. In particular the 
size of the visible (causal) horizon today is ~ 3000 Mpc. 

• The fossil Cosmic Microwave Background (CMB) radiation: The Universe is immersed in a bath of 
thermal photons at a temperature To — 2.73K with an almost perfect blackbody distribution. This 
distribution and small anisotropies of order AT/To ^ 10^^ were measured in 1992 by the COBE 
satellite and their detection represents a triumph for the standard hot big bang model The small 
temperature anisotropies, predicted by cosmological models, provide the clue to the origin of galaxy 
formation and large scale structure and is an important confirmation of theories of early Universe 
cosmology. 

• The abundance of light elements: Observations of the abundance of elements in low metallicity regions 
reveals that about 76% of matter is in the form of Hydrogen, about 24% (by mass) in "^He and very 
small abundances oi^He (~ 10~^), deuterium (D) (~ 10"^) and Li (~ 10-^°) all relative to Hydrogen. 
These elements were formed during the first three minutes after the Big Bang, while heavier elements 
are produced in the interior of stars and astrophysical processes such as supernovae explosions. 



The main building blocks for a theory of the standard Hot Big Bang are: 

• Gravity: Classical general relativity provides a good description of the geometry of space time for 
distances I > Ipi ^ 10~'^^cm or time scales t > tpi ^ 10~'*'^s, or equivalently energy scales smaller than 
the Planck scale Mpi ~ lO^^Gev. We have to wait for a consistent quantum theory of gravity unified 
with matter to explain phenomena on shorter space-time scales or larger energy scales. 

Homogeneity and isotropy lead to the Robertson- Walker metric 



where t is the comoving time. The constant k determines the spatial curvature and can be set to be 
either ±1 or by redefining the scale of coordinates. For k — 1,0,-1 the Universe is closed, flat or 
open respectively. The scale factor a{t) relates physical and comoving scales 



2.2 The building blocks: 




(2.3) 




a{t)l. 



corn 



(2.4) 



The Friedmann equation determines the evolution of the scale factor from the energy density 




(2.5) 



A spatially flat Universe has the critical density 
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Pc - ^^liH^ (2-6) 

and it is customary to introduce the ratio of the density of any component (radiation, matter etc) to 
the critical density as 

17 = — (2.7) 

Pc 

The energy momentum tensor is assumed to have the fluid form leading to the (fluid) conservation of 
energy equation 

p + 3(p + p)-=0 (2.8) 
a 



where p, p are the energy density and pressure respectively. The two equations ( p.5y2.8| ) can be combined 
to yield the acceleration of the scale factor, 



^{p + 3p) (2.9) 



which will prove useful later. In order to provide a close set of equations we must append an "equation 
of state" p = p{p) which is typically written in the form 



p^wip)p (2.10) 

Current observations favor a spatially flat Universe fc = which is consistent with predictions from 
inflationary scenarios. For a spatially flat Universe k = and for w — constant we obtain the following 
important cases: 



w = : Matter domination =^ p cx a ^ ; a{t) cx (2-11) 

w — — : Radiation domination ^ p oc a^^ ; a{t) oc i ^ (2-12) 
3 



w = — 1 : De Sitter expansion => p = constant ; a{t) oc e^^ : H = . (2-13) 

Furthermore we see that for accelerated expansion it must be the case that w < —1/3. 

• The Standard Model of Particle Physics: the current standard model of particle physics, experi- 
mentally tested with remarkable precision describes the theory of strong (QCD), weak and electromag- 
netic interactions (EW) as a gauge theory based on the group SU{3)c ^ SU{2) (g) U{1). The particle 
content is: 3 generations of quarks and leptons: 




vector Bosons: 8 gluons (massless) , Zq (mediate neutral currents), (mediate charge currents) with 
masses of order 80, 90 Gev respectively and the photon (massless) and scalar Higgs bosons, although 
the experimental evidence for the Higgs boson is still inconclusive. 
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Current theoretical ideas propose that the strong, weak and electromagnetic interactions are unified in 
a grand unified quantum field theory (GUT), perhaps with supersymmetry as the underlying fundamental 
symmetry and a unification scale Mqut ^ lO^^Gev. Furthermore the ultimate scale at which Gravity is 
eventually unified with the rest of particle physics is the Planck scale Mpi ^ lO^^Gev. Although there are 
several proposals for the total unification of forces some of which invoke strings, M-theory, extra dimensions 
and a variety of novel and fascinating new concepts, the experimental confirmation of any of these ideas 
will not be available soon. However, the physics of the "standard" model of the strong and electroweak 
interactions that describes phenomena at energy scales below ~ 100 Gev is on solid experimental footing. 

The connection between the standard model of particle physics and early Universe cosmology is through 
Einstein's equations that couple the space-time geometry to the matter-energy content. We argued above 
that at energy scales well below the Planck scale gravity can be studied classically. However, the standard 
model of particle physics is a quantum field theory, thus the question arises: how to treat classical space-time 
but with sources that are quantum fields. The answer to this question is: semiclassically, Einstein's equations 
(without the cosmological constant) are interpreted as 



g,M">^^M>'__^gM">^=__(yM<>) (2.14) 
2 Mpi 

with r^'' is the operator energy momentum tensor and the expectation value is taken in a given quantum 
state (or density matrix). A state that is compatible with homogeneity and isotropy must be translational 
(and rotational) invariant, and the expectation value of the energy momentum tensor operator must have 
the fluid form (T^) = diag[p, — p, — p, —p]. 

Therefore, through this identification the standard model of particle physics provides the sources for 
Einstein's equations. All of the elements are now in place to understand the evolution of the early Universe 
from the fundamental standard model. Einstein's equations determine the evolution of the scale factor, the 
standard model provides the energy momentum tensor. 



2.3 Energy scales: 

While a detailed description of early Universe cosmology is available in many excellent books a broad- 

brush picture of the main cosmological epochs can be obtained by focusing on the energy scales of particle, 
and atomic physics. 

Total Unification: Gravitational, strong and electroweak interactions are conjectured to become uni- 
fied and described by a single quantum theory at the Planck scale ~ lO^^Gev. While there are currently 
many proposals that seek to provide such fundamental description, these are still fairly speculative and no 
experimental confirmation is yet available. 

Grand Unification: Strong and electroweak interactions (perhaps with supersymmetry) are conjectured 
to become unified at an energy scale ~ 10^^ Gev corresponding to a temperature T ~ IQ'^^K. There are 
very compelling theoretical reasons (such as the joining of the running coupling constants) that lead to this 
conjecture, but there is as yet no experimental evidence in favor of these ideas. 

Electroweak : Weak and electromagnetic interactions become unified in the electroweak theory based 
on the gauge group SU{2) (g) U{1)y- The weak interactions become short ranged after a symmetry breaking 
phase transition SU{2) (g) U{1)y U{l)em, at an energy scale of the order of the mass of the Zo,W^ 
vector bosons, E ^ lOOGev corresponding to a temperature Tew 100 Gev. At temperatures T > Tew 
the symmetry is restored and all vector bosons are (almost) massless (save for plasma effects that induce 
screening masses) . For T < Tew the vector bosons that mediate the weak interactions (neutral and charged 
currents) , Zq acquire masses while the photon remains massless. Thus Tew determines the temperature 
scale of the electroweak phase transition in the early Universe and is the earliest phase transition that is 
predicted by the standard model of particle physics. 
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QCD: The strong interactions have a typical energy scale Kqcd ~ 200Mev at which the coupling 
constant becomes of order one. This energy scale corresponds to a temperature scale Tqcd ^ lO^^K. 
QCD is an asymptotically free theory, the coupling between quarks and gluons becomes smaller at large 
energies, but it diverges at the scale Aqcd- For energy scales below Aqcd QCD is a strongly interacting 
theory.This phenomenon is interpreted in terms of a phase transition at an energy scale Aqcd ~ 200 Mev 
or Tqcd ^ 10^^ K. For T > Tqcd the relevant degrees of freedom are weakly interacting quarks and gluons, 
while below are hadrons. This is the quark-hadron or confinement-deconfinement phase transition. At about 
the same temperature scale QCD has another phase transition that results in chiral symmetry breaking (for 
more details see section 4). The QCD phase transition(s) are the last phase transition predicted by the 
standard model of particle physics. The high temperature phase above Tqcd, with almost free quarks and 
gluons (because the coupling is small by asymptotic freedom) is a quark-gluon plasma or QGP for short. 
Current experimental programs at CERN (SPS-LHC) and Brookhaven National Laboratory(AGS-RHIC) 
are studying the QCD phase transition via ultrarelativistic heavy ion collisions and a systematic analysis of 
the data gathered at SPS-CERN during the last decade has given an optimistic perspective of the existence 
of the QGpH, H (see section 4 below). 

Nuclear Physics: Low energy scales that are relevant for cosmology are determined by the binding 
energy of light elements, in particular deuterium, whose binding energy is ~ 2 Mev corresponding to a 
temperature T^s ~ 10"'^'' K. This is the energy scale that determines the onset of primordial nucleosynthesis, 
as described below. 

Atomic Physics: Another very important low energy scale relevant for cosmology corresponds to the 
binding energy of hydrogen ^ 10 eV. This is the energy scale at which free protons and electrons combine 
into neutral hydrogen. As it will be described below the relevant scale is more like < 0.3 eV corresponding 
to a temperature T ~ 3 x lO'^ K. 

Based on the ingredients described above, a very detailed picture of the thermal history of the Universe 
emerges ||-§: during the first - 10000 years after the Big Bang the Universe was radiation dominated 
expanding and cooling (almost) adiabatically. As a consequence the entropy S oc V{t)T^{t) cx Vq [a{t)T{t)f' 
is constant, implying T — To/a(t). Radiation domination, in turn, results in that a{t) cx and a detailed 
analysis reveals that 



m ^ 4^ ^ 1^ (2.15) 

t2[sec) t^{sec) 

3 Phase Transitions and their aftermath 
3.1 GUT'S and inflation 

Current theoretical ideas of theories beyond the standard model suggest that there could have been a phase 
transition at the GUT scale Tqut ^ 10^^ Gev ^ lO'^^ K. This energy scale is also usually associated with 
an important cosmological stage: inflation, during which the acceleration of the scale factor is positive which 



implies via eqn. (2^) that p < —p/i- Inflation plays a very important role in early Universe cosmology and 
current observations of the power spectrum of temperature anisotropies seem to confirm the robust features 
of the inflationary proposal ||l[ . While I will not attempt to review all features of inflationary cosmology, for 
which the reader is referred to the literature I summarize some of the most important concepts so as 

to make contact with the observable consequences of phase transitions later. 

As is mentioned above, inflation corresponds to an epoch of accelerated expansion, i.e, with a{t) > 



which from eqn. (2.9) requires that p < — jo/3. Within particle physics models this is achieved by considering 
the energy momentum tensor of a scalar field, which in principle is one of the fields in the GUT. For a scalar 
field the energy density and the pressure are given by [^-0 
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p 
p 



(3.1) 
(3.2) 



where 1^(1^9) is the scalar potential. 

Inflation results from the generalized slow- roll condition §-0 



which results in that p = —p ~ V{(p) ~ constant leading to the De Sitter solution cqn.( ^.13 ) for the scale 



factor. This situation can be achieved via a variety of inflationary scenarios (old, new, chaotic, hybrid 
etc,) With the purpose of establishing contact with observational consequences of phase transitions I 

focus the discussion on either a first order (old inflation) or second order (new inflation) phase transition 
at the GUT scale. 

In these situations the expectation value of the scalar field is at a false vacuum extremum of the potential 
with a very slow time evolution. The expansion of the Universe red-shifts the inhomogeneities of the field 



and the slow-roll condition (3.2) is fulfilled. In these scenarios the expectation value of the scalar field is 
nearly constant during the infiationary stage leading to a De Sitter expansion of the scale factor as in eqn. 
( ^.13D with 



H^J'-^ (3.4) 

Inflation terminates when the field "rolls down" to the minimum of potential, where further oscillations 
result in particle production, reheating and a transition to a standard hot big bang, radiation dominated 
era While the evolution of the scale factor along with the dynamics of the scalar field (inflaton) 

were typically studied using the classical equations of motion for the scalar field, more recently a consistent 
description of the quantum dynamics has been provided ^ [T^. The non-equilibrium dynamics of cos- 
mological phase transitions requires a non-perturbative framework that treats self-consistently the dynamics 
of the matter field and the evolution of the metric. This framework leads to a detailed understanding of 
the classicalization of quantum fluctuations as well as to a microscopic justification for inflation in a full 
quantum field theory. We refer the reader to ref. |[ for more details which fall outside the main scope 
of these lectures. 

An important aspect of inflationary dynamics and one that is very relevant to the discussion of observables 
associated with phase transitions stems from the positive acceleration of the scale factor. A simple calculation 
shows that 



^>0=^^>^ (3.5) 
a{t) a{t) dnit) 

Therefore during a period of accelerated expansion or inflation, the scale factor grows faster than the 
Hubble radius. In particular for De Sitter inflation while the Hubble radius dn = ^/H is constant, the 
scale factor grows exponentially a{t) = age^*. This feature of accelerated expansion is indeed remarkable. 
Consider a perturbation of physical wavelength Xp^ysit) = Aca(t) where Ac is the comoving wavelength. 
When Xphys{t) < da causal microphysical processes can affect this perturbation, but when this wavelength 
"crosses the horizon" i.e, when Xphys{t) > dn no causal process can affect the perturbation. Thus physical 
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wavelengths inside the Hubble radius (or loosely speaking the causal horizon) are causally connected and 
influenced by microphysical processes. When these wavelengths cross the horizon they stretch superluminally 
and "decouple" from causal processes, hence their evolution is not affected by microphysical processes. 

We note that for matter or radiation dominated cosmologies a{t) cx t^^"^ , t^/^ respectively, for which the 
Hubble radius dnit) = 3f/2 , 2t respectively, grows faster than the scale factor. Therefore wavelengths that 
at some time are outside the Hubble radius, eventually "cross the horizon" back inside the Hubble radius. 

Since inflation is followed by the standard hot big bang cosmology with a radiation and matter dominated 
eras, physical wavelengths that are inside the Hubble radius during inflation, and cross outside during 
inflation will re-enter during the radiation or matter dominated eras. This situation is depicted in figure (|l|) 
below. 



\ (Today's Hubble radius) 




3eq ^0 Ln(a) 



Figure 1: Logarithm of physical scales vs. logarithm of the scale factor. The causal horizon dn is shown for 
the inflationary (Dc Sitter), radiation dominated and matter dominated stages. The physical wavelengths 
i^phys) for today's Hubble radius, the surface of last scattering (LSS) and a galactic scale are shown. Qeq , ao 
refer to the scale factor at matter-radiation equality and today, respectively. 



The wavelengths that are of cosmological relevance today say with 1 Mpc < Xphys{ta) S 3000 Mpc crossed 
the horizon during the last ~ 10 e — folds of inflation|^. 

3.1.1 Density perturbations and the signature of a phase transition 

This important feature of accelerated expansion and inflation, i.e, that physical wavelengths within the 
causal horizon during inflation cross-out and re-enter during radiation or matter domination, provides the 
mechanism for "seeding" temperature anisotropies. Inflation "seeds" temperature inhomogeneities in the 
cosmic microwave background and the matter density inhomogeneities that lead to large scale structure 
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formation from primordial quantum fluctuations whose wavelengths were inside the Hubble radius during 
inflation§-§. 

To see this more clearly, consider small amplitude perturbations of the homogeneous scalar field that 
drive inflation 



Lp{x,t)^ip{t)+5ip{x,t) (3.6) 

This perturbation induces a small perturbation in the energy density, which to linear order in Sip(x, t) 
and taking the spatial Fourier transform is (3|-0| 

h^^-y^ 5n (3-7) 

The power spectrum of density perturbations is obtained from the quantum average 

(l%P> « (3.8) 
An important and robust result from inflation is that 

oc fc"=-i ; with n, - 1 (3.9) 

This power law spectrum with rig ^ 1 is (almost) scale invariant and referred to as Harrison-Zeldovich. 
As it will be discussed in the section on the CMB below, this power spectrum is measurable through the 
temperature anisotropics at the surface of last scattering |^ |j . 

The main point of this discussion is that inflation generates temperature anisotropics and density per- 
turbations from quantum fluctuations whose wavelengths cross the horizon during inflation and re-enter just 
before recombination ^ ^ (see fig. p. 

Phase transitions during infiation modify the power spectrum of the quantum fiuctuations. After the 
wavelength of these fluctuations crosses the horizon during inflation, they evolve acausally carrying with 
them the information of the phase transition without being affected by microphysical processes. When 
the wavelengths of these fluctuations re-enter the horizon near recombination the information of the phase 
transition is imprinted in the temperature anisotropics in the CMB through departures from scale invariance 



in the power spectrum (3.8), i.e, ^ 1 as a consequence of the phase transition. 



To illustrate these ideas in a simple setting, let us consider a simple scalar field theory with the potential 



= (3.10) 

This is the typical potential that leads to (a second order) phase transition. The strength of the self- 
couphng A is constrained by the amplitude of the temperature anisotropics to be [|| A ~ 10^^^. During the 
phase transition, the expectation value of the scalar field ^ 0. Requiring that the energy density is given 
by the GUT scale 

p-^-(10^^ Gev)4 

during the epoch of the phase transition determines that m ^ 10^^ Gev. The Hubble constant is then given 
during this stage by 
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Therefore during the phase transition when (p ^ the expansion is of De Sitter type and small amplitude 
fluctuations of the scalar field (in terms of the spatial Fourier transform) Sipj: obey the linearized equation 



6^At) = 0. 



(3.12) 



We note that the sign of the mass squared term is negative as required to describe a symmetry breaking 
phase transition. Setting H — Q , a — \, i.e, in Minkowsky space time we recognize that long- wavelength 
fluctuations with fc^ < vn? grow exponentially. These are the spinodal instabilities, a hallmark of the process 



of phase separation during the phase transition jTo[ |T2j. During the stage of De Sitter inflation, eqn.( 3.1^ ) 
has the following solution |9) 



S<ff:{t) = e" 



-Ht 



BkJ-p ( — e 



-Ht 



(3.13) 



the coefficients Ak,Bk are determined by initial conditions]^. When the physical wavevector crosses the 
horizon, i.e, when k e^^' ^ H the solution grows exponentially and is given by 



Sipf:it) oc e('^-i)^* (3.14) 

In the limit H « m this solution displays the spinodal long- wavelength instabilities cx e™* associated 
with the phase transition in Minkowski space-time jl^, |l^. Thus the exponential growth of superhorizon 



fluctuations given by (3.14) is a direct consequence of the spinodal instabilities associated with the phase 



transition during the De Sitter inflationary epoch. A detailed calculation of the power spectrum of density 
fluctuations |ll[ yields 



cx fc-2(^-i) ^ns = l~2{v-^) (3.15) 

since v > 3/2 we see that there is more power at long wavelengths. This is obviously a consequence of 
the spinodal instabilities associated with the phase transition. Thus we arrive at an important result that 
provides an observational signal of the phase transition: the power spectrum of density perturbations is 
tilted to the red, i.e, the index is smaller than one. This is because the spinodal instabilities associated 
with the phase transition imply larger amplitudes for longer wavelengths, thus enhancing the infrared. After 
these perturbations re-enter the horizon close to recombination, this power spectrum is imprinted in the 
temperature anisotropics of the CMB (see the discussion on the CMB below) which are measured. The four 
year COBE-DMR Sky Map[|l3j gives Ug ^ 1.2 ±0.3, thus a (second order) phase transition at the GUT scale 
is consistent with this powe r spectrum if H/m > 3 which, in turn, is consistent with GUT scale inflation as 
can be seen from eqn.( |3.11 ). 

This simple example thus provides a good idea of the potential observational signatures of phase transi- 
tions at the GUT scale associated with an inflationary era: a power spectrum with a red-tilt results if the 
phase transition is triggered by spinodal instabilities during the inflationary epoch. 
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Another potential observable from an inflationary phase transition was proposed recently p4| . These 
authors pointed out that the breaking of scale invariance in the primordial power spectrum by an inflation- 
ary phase transition could lead to a step-like spectral feature, which seems to be compatible with recent 
measurements of the CMB and large scale surveys (APM). 

Thus observations of the temperature anisotropics in the CMB could provide information on inflationary 
phase transitions. Furthermore, a "reconstruction" program seeks to extract some aspects of the scalar 
potential from the CMB temperature anisotropics [^5j, therefore providing a further link between phase 
transitions during the inflationary era and the temperature anisotropics of the CMB. 



3.2 The Electroweak Phase Transition: Baryogenesis? 

The EWPT is the first symmetry breaking phase transition predicted by the standard model of particle 
physics. It occured at a temperature Tew ~ 100 Gev ^ lO^^K at about t ~ 10~^^ sees after the Big 
Bang when the Hubble radius was djj ^ 10~^cm. The symmetry breaking pattern is SU{2) (g) U{1)y 
U{l)em- Probably the most tantalizing observable from the electroweak phase transition could be the baryon 
asymmetry. There is an asymmetry between particles and antiparticles in the observed Universe, large regions 
of antimatter would result in particle annihilations leading to a diffuse 7-ray background and distortions of 
the CMB. None of which is observed, leading to the conclusion that the Universe is made up of particles 
up to the Hubble radius iQ. Furthermore Big Bang nucleosynthesis provides accurate predictions for the 
abundance of light elements up to ^Li in terms of the ratio of the baryon density to the photon density 

rj=!^^ (3.16) 
Observations of the abundance of light elements constrain this parameter to be in the range |^ 

4 X 10"^° < 77 < 7 X 10"^° (3.17) 

The important question is: what is the origin of the baryon asymmetry (baryogenesis)? (for a recent 
review on baryogenesis see JTs]). 

The necessary conditions required for baryogenesis were identified originally by Sakharov [ pTf 

• C and CP violation 

• Baryon number violation 

• Departure from thermal equilibrium 

The EW theory violates maximally P (parity) and also C (charge conjugation). CP violation is observed 
in the K^,K^ system (and currently also in the B^,B'^ system) and is experimentally determined by the 
parameters e ^ 2.3 x 10"'^, e' ~ 5.2 x 10"^ that measure indirect (through mixing of eigenstates) and direct 
CP violation. CP violation in the standard model is a consequence of the phase in the CKM mass matrix 
for three generations. 

Baryon number violation in the EW theory is a consequence of non-perturbative gauge (and Higgs) field 
configurations that interpolate between topologically different vacua. A transition between two adjacent 



vacua leads to one unit of baryon number violation per family (for details see |16 ). While at zero temperature 

27r 

this is a vacuum tunneling process which is suppressed by the barrier penetration factor e °™ with ^ 
1/30, at finite temperature the transition is overbarrier, i.e, thermal activation. The overbarrier transitions 
are unsuppressed at temperatures of the order of Tew [0 • 

The requirements of baryon number violation and CP violation are clear: to generate any asymmetry the 
theory must allow for processes that produce the quantum number in question (B) and to also distinguish 
particle from antiparticle (CP). If the relevant processes are in equilibrium then the rate for particles and that 
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of antiparticles are related by detailed balance and since CPT invariance implies that the masses of particles 
are equal to those of antiparticles, the distribution functions for both are identical and no asymmetry can 
result g 0. 

The expansion of the Universe alone cannot provide the non-equilibrium ingredient. To see this let us 
compare the reaction rate of weak interactions (the electromagnetic interactions are much faster than the 
weak interactions) to the expansion rate of the Universe. Since the EWPT takes place during a radiation 
dominated era with p oc the expansion rate is given by 

i/oc|!- (3.18) 

The reaction rate is given by F = na^ with n the particle density and ct„ the cross section. The typical 
weak interaction cross section is ~ G'^E^ with Gp ~ 10^^/Gev^, at temperature T the typical energy is 
E ^ T and a particle density n ^ . Thus the weak interaction rate T = naw is of order 



therefore the ratio 



r 

H 



GlT'^ (3.19) 



GlMpiT^ ^ [lO^T{Gev)]^ (3.20) 



determines whether the weak interactions are in local thermal equilibrium (LTE). At T ^ Tew the ratio is 
T/H ^ 10^^, therefore weak interaction processes are in LTE (and obviously electromagnetic processes will 
be more so). 

Instead departures from equilibrium results from a strong first order phase transition that occurs via the 
nucleation of bubbles of the true phase in a background of the metastable phase ||l^ . 

Detailed numerical simulations of the EWPTjlSf reveal that it is first order if the Higgs mass Mh < 80Gev 
but above this value it is a smooth crossover. Recently, however, a summary of all the data collected by the 
four detectors at LEP up to energy 202 Gev, seem to indicate that the Higgs has a mass Mh ^ 115Gev[p^. 
Thus it is very likely that the (minimal) standard model cannot accomodate a strong first order phase 
transition that will produce the non-equilibrium conditions sufficient for baryogenesis. Furthermore, it 
is now clear that the CP violation in the (minimal) standard model, encoded in the parameters e, e' (or 
alternative in the phase of the CKM matrix) is too small to explain the observed baryon asymmetry. 

Thus, the current theoretical understanding seems to suggest that the baryon asymmetry cannot be 
explained by the minimal standard model although the Minimal Supersymmetric Standard Model (MSSM) 
may accomodate all of the necessary ingredients [^6). Hence the cosmological consequences of the EWPT 
are not very clear. Furthermore, while accelerator experiments (Tevatron and LHC) will probably lead to 
an assessment of the Higgs mass and CP violating parameters, a direct study of the EWPT will not be 
feasible soon. Indeed in order to study it an energy density of order pEW '^biv ™ust be achieved, but 
with Tew ^ 100 Gev this energy density is pew lO^^Pn where p„ ^ 0.15 Gev/fm'^ is the energy density 
of nuclear matter!!. 

3.3 The QCD phase transition(s) 

The next and last phase transition predicted by the standard model of particle physics takes place at the 
QCD scale T ~ 100 — 200Mev. It occurred when the Universe was about t ^ 10~^ — 10~^ seconds old and the 
Hubble radius was dn ^ lOKm. QCD is an asymptotically free theory, that is the gauge coupling constant 
between quarks and gluons varies with energy through "vacuum polarization" effects as 
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as{E)^ f ; AQCD-- 200Mev ; a, = ^ (3.21) 



(11-1^/) In 



4tt 



with Nf being the number of flavors and g is the quark-giuon couphng. For E 3> ^qcd quarks and gluons are 
weakly coupled, but for E ^ Aqcd the couphng becomes strong, diverging at i? = ^qcd- This divergence 
is interpreted as a transition between a state of almost free quarks and gluons to a state in which quarks 
and gluons are confined inside hadrons. This is the confinement- deconfinement transition and the QCD scale 
Aqc d determines the range of temperature for this transition. The high temperature phase in which quarks 
and gluons are almost free is referred to as the Quark Gluon Plasma or QGP for short. 

Because the up and down quark masses are so much smaller than Aqcd {mu 5 Mev; ~ 10 Mev) the 
low energy sector of QCD has an SU{2)ji® SU{2)l symmetry in the limit of massless up and down quarks. 
This chiral symmetry corresponds to independent chiral rotations of the right and left handed components 
of the quark fields. However the ground state of QCD spontaneously breaks this chiral symmetry down to 
S'[/(2)l+_r leading to three massless Goldstone bosons. The small (and different) masses of the up and down 
quark provide a small explicit symmetry breaking term in the Lagrangian and as a consequence the (would be) 
Goldstone bosons acquire a mass. This triplet of almost Goldstone bosons are the neutral and charged pions, 
which are the lightest pseudoscalar particles with m^r ~ 140 Mev. Lattice simulations |2^ suggest that the 
confinement-deconfinement and chiral phase transitions occur at about the same temperature T ~ 160 Mev 
(for a review see [^ll ). At this temperature only the lightest quark flavors influence the thermodynamics, 
the up, down and strange quarks with ^ 150 Mev. The chiral phase transition associated with the two 
lightest quark flavors is second order for massless quarks but either slightly first order or a crossover when 
the mass of the up and down quarks is accounted for, and is in the universality class of the 0(4) Heisenberg 
ferromagnet. Including the strange quark the situation is more complicated and depends on whether the 
strange quark can be considered heavy or light, the transition becomes first order if the strange quark is 



light 1 20 . The complication arises because rus ^ Aqcd and is therefore neither light nor heavy on the QCD 
scale. 

The confinement-deconfinement transition for three flavors is likely to be of first order |2|] and 
therefore occurs via the nucleation of hadronic bubbles in the background of a QGP. An estimate of whether 
the expansion of the Universe at the QCD phase transition results in non-equilibrium effects can be obtained 
in the same manner as in the case of the EW transition, i.e, by estimating the ratio T / H with F a typical 
strong interaction reaction rate. A typical cross section is u '--^ 0.1 — 1 fm^ and the number of particles at 
Tqcd ^ 200 Mev is n - T^cd ^ Vfm^ this leads to F - 1/fm/c - 10^^ sec~^ and since H - lO*' sec-\ 
therefore V /H ~ 10^^ and again the QGP is in LTE (local thermal equilibrium). However, just as in the case 
of the EWPT, non-equilibrium effects may arise from a first order phase transition through supercooling 
and the nucleation of bubbles. In this case, hadronic bubbles nucleate in the host of the QGP releasing 
latent heat. After a short period of reheating, the transformation between the QGP and the hadronic phases 
occurs in LTE until all of the QGP transforms into hadrons. Most of the hadrons decay on short time scales, 
for example charged pions decay on time scales ~ 10~^ sees, and neutral pions on much shorter time scale 
10~^^ sees. Therefore most hadrons have decayed after a few /^secs after the QCD PT but for neutrons, with 
a lifetime ~ 900 sees and protons (with a lower limit for the lifetime of ^ 10^^ years). 

Since the QCD PT is the last phase transition, there is an important effort to understand its potential 
cosmological consequences. While the details are complex, not completely understood and cannot be given 
in this short review, the following important consequences are possible (for a recent review see also 22, 23| ) 

• Baryon inhomogeneities affect nucleosynthesis: ||2j] The nucleation of hadronic bubbles in a host 
of QGP can lead to inhomogeneities in the baryon number density, with the scale for inhomogeneity 
determined by the typical distance between bubbles. These inhomogeneities may result in inhomo- 
geneous neutron to proton ratios which in turn can lead to inhomogeneous nucleosynthesis (see the 
section on nucleosynthesis below) and modify the abundance of light elements. For the inhomogeneities 
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produced by bubble nucleation to modify the neutron to proton ratio for nucleosynthesis it must be 
that the typical separation between nucleating bubbles must be much larger that the proton and neu- 
tron diffusion lengths. Proton and neutrons have different diffusion lengths because of the Coulomb 
scattering of protons. An important criterion for inhomogeneous nucleosynthesis is that the mean 
scale for baryon inhomogeneity at the time of the QCD PT, i.e, the mean distance between nucleating 
bubbles must be of order ~ 1 meter . This distance depends on the latent heat released during the 
(first order) phase transition, the free energy difference between the QGP and the hadron gas phases 
and the surface tension for the nucleating bubble. More recent analysis seems to point out that 
the distance between bubbles produced during homogeneous nucleation is dnuc 1cm. Larger values 
of dnuci can arise from inhomogeneous nucleation seeded by "impurities" [E5|. These inhomogeneities 
can be produced by primordial black holes (see below) as proposed by ||25| or by primordial density 
fluctuations as suggested by |p6|| . 

However as recognized in the literature |24, 2^, ^ much more needs to be understood about the phase 



structure of QCD before a quantitatively reliable statement can be made. Of particular importance 
are details of the equation of state (Eos), surface tension, which determines the size of the nucleat- 
ing bubbles, and also the (in medium) mass of the strange quark. For a more recent estimate on 



inhomogeneous nucleosynthesis arising from a first order QCD PT see 27, p3| 



Primordial Black Holes: An important aspect of a first order phase transition (with only one 

globally conserved quantity, such as baryon number) is that during the coexistence or mixed phase of 
QGP and hadron gas, the Gibbs construction determines that the phase transition occurs at constant 
pressure (as well as temperature). If the pressure is constant, then the (adiabatic) speed of sound 
vanishes and there are no restoring pressure waves. This in turn means that there is no restoring force 
to counterbalance the gravitational collapse and long- wavelength density perturbations will grow under 
self-gravity causing the gravitational collapse of the mass contained in a volume of radius given by this 
wavelength. This observation leads to the possibility of formation of primordial black holes p8|(for 
a discussion of gravitational collapse and Jeans instability see section 3.6). An important question is 
what is the mass of the black hole?, to obtain a qualitative estimate we can calculate the mass in the 
Hubble radius at the time of the QCD PT: Mh = 47rp /3 with dn ~ lOKm the Hubble radius at the 
time of the PT. Using that = l/d% 8np /iMl^ therefore Mh = MlidHl2 ~ lO^Gev ~ IMq. 
Thus this order-of-magnitude estimate suggests that primordial black holes with mass up to ~ IM© 
may form during the QCD PT |^. It is important to mention at this stage, that the Hubble radius 
at the QCD PT ~ 10'' m extrapolated to today is ~ 0.2 pc which is much smaller than the Hubble 
radius today ^ 3000 Mpc. Therefore a large number of primordial black holes produced during the 
QCD PT would be present today. 

A more detailed analysis of the possibility of primordial black hole formation and the ensuing density 
fluctuations has been provided recently H^. The results of this reference, based on approximate or 
lattice EoS seem to lead to a less likely scenario for solar mass primordial black hole formation, but 
to an enhancement in the clumping of cold dark matter (CDM). It is fair to say, however, that there 
are still large uncertainties in the EoS and the relevant quantities that enter in the calculation. The 
experimental programs that seek to study the QGP in heavy ion colliders that are described below will 
lead to a more reliable understanding of the relevant aspects of the PT, EoS and QCD parameters. 
For more detailed aspects of primordial black holes, see pof . 

Strange Quark Nuggets: Witten suggested that during a first order QCD PT in the early 
Universe a large amount of strange matter could be produced He also argued that strange 

matter could be absolutely stable with an energy per baryon which is less than the maximum binding 
energy (of Fe) 930Mev. While this scenario has been criticized (see Applegate and Hogan in ) , the 
uncertainties in the knowledge of the QCD EoS, the (in medium) strange quark mass and other relevant 
QCD parameters leave enough room for the possibility of formation of strange quark "nuggets" . 
Some of the consequences of the formation of strange quark nuggets during the (first order) QCD PT 
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had been investigated in These strange nuggets can be part of the cold dark matter, and if 

they do not evaporate during ~ 1 sec after the QCD PT, their presence can affect nucleosynthesis by 
modifying the neutron to proton ratio, thus modifying the "^He (and heavier elements) abundance 

The different possible signatures of the QCD PT in the early Universe as described by the scenarios above 
rely on the particular details of the phase transition as well as thermodynamic parameters that cannot be 
calculated in perturbation theory. Lattice gauge theory provides a non-perturbative approach to studying 
many of these aspects and although progress has been made in the field the complications associated with 
light quarks beyond the quenched approximation still need further understanding, ft is clear that probably 
the best assessment of the QCD phase diagram would be obtained from experiments that can probe the hot 
and dense phase(s) of QCD. 

The heavy ion programs at the CERN-SPS and the BNL-AGS, and the current program at BNL-RHIC 
have as major goals to reveal the new state of matter described by the QGP. Almost two decades of experi- 
ments at CERN and BNL have provided a wealth of data that together reveal that this new state of matter, 
predicted by QCD may have been formed in ultrarelativistic heavy ion collisions. The recent announcements 
from CERN jsj, ^ yield convincing arguments that there is already experimental evidence for the QGP. 
Thus there is the very tantalizing possibility that the next generation of ultrarelativistic heavy ion collisions 
at RHIC and eventually the ALICE (A Large Ion Collider Experiment) program at LHC will provide a more 
detailed understanding of the QCD PT, and from that data we can learn about the observable consequences 
of the same PT in the early Universe. This program, along with a brief summary of the recent results will 
be described in more detail in section 4. 

3.4 Nucleosynthesis 

The next stage down the ladder of energy scales corresponds to ~ 1 Mev. This is the scale of binding energy 
of deuterium and determines the onset of primordial nucleosynthesis. As mentioned in the previous section, 
after the QCD confinement-deconfinement phase transition, the thermodynamics is described by a hadron 
gas. After about a /isec after this PT most hadrons decay and only neutrons and protons remain. Therefore 



the Universe can be characterized by a plasma of n,p, e*, ly's, 7's. Equation (3.20) for the ratio of a typical 
weak interaction reaction rate and the expansion rate of the Universe, determines that at a temperature 
T ~ 1 Mev the weak interactions freeze-out (decoup le), i. e, the rate for weak interaction processes is smaller 



than that of the expansion of the Universe. Eqn. (3.20) also means that the mean free path for the weak 



interactions is larger than the Hubble radius for T < 1 Mev. A more detailed calculation ||36|] shows that 
the freeze-out or decoupling temperature is actually Tf = 0.8 Mev. For T ^ Tf neutrons, protons, electrons 
and neutrinos are in nuclear statistical equilibrium via the weak interaction reactions 

n ^ p + + De '■, Ve + n ^ p + ; e'^ + n ^ p + (3.22) 

Since neutrons and protons are non-relativistic near the freeze-out temperature, their phase space dis- 
tribution is determined by the usual Boltzmann distribution function and at any given temperature much 
smaller than their masses the neutron to proton ratio is given by 

Nr, Am 1.29 

— = e"— = e"''''^ (3.23) 
with Am the neutron-proton mass difference. The neutron to proton ratio at T ^ T/ is given by 

^ ^g-1.29/0.8_^/g (3.24) 
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After the Universe has cooled to the freeze-out temperature Tf this ratio remains fixed but for the decay 
of the neutron with a hfetime ^ 900 sees. 
Nucleosynthesis begins via the reaction 

n + p^ d + -f (3.25) 
and continues through the following two possible paths towards the conversion of all neutrons into "^He 



d + d'^'-^He + n ; or d + d<^^H +p 
^He + d <-> ^He +p ; or ^H + d^ ^He + n (3.26) 

The binding energy of deuterium is 2.22 Mev and a priori one would think that once the temperature 
falls below this value the reaction ( 3.25| ) will mainly go to the right. However, because the number of photons 



is so much larger than the number of baryons since ng jn^ cx r] with r] given by eqn. ( ^.17 ) there are too 



many photons in the high energy tail of the blackbody distribution to photodisintegrate any deuterium formed 
until the temperature actually falls to T ^ 0.1 Mev at a time of about 2 minutes P|. This is the "deuterium 



bottleneck": nucleosynthesis does not begin in earnest until the first step (3.25) happens. The formation 
of deuterium is hindered by the fact that the Saha equation for the equilibrium abundance for a species 
of atomic number A, contains a factor cx rj^^^ From the time of freeze-out when the temperature was 
T Tf until deuterium is formed, the neutrons have decayed via the weak interactions (with a lifetime 900 
sees) to a neutron to proton ratio Nn/Np ~ 1/7 [p| As soon as deuterium is formed the second stage(s) 



( f3.26 ) occur fairly fast ||, ^ and all neutrons are used up in the formation of ^He. The neutron to proton 



ratio determines the abundance of ^iJe: each atom of '^He has two protons and two neutrons, but for each 
two neutrons 14 protons are needed because of the ratio 1/7. Altogether there are 16 baryons, out of which 
one Helium atom was formed, so the abundance of Helium is 25%!. A network of nuclear reactions leads to 
the formation of light elements up to ^Li |3^. However because there is no nucleus with A — 5 and no stable 
one with A — 8 with unstable nuclei in between i^Li^'^ Be) only small quantities of light elements other than 
^He are produced. An important ingredient in these calculations, besides the nuclear cross sections that 
determine the nuclear reaction rates that have been measured in low energy experiments is the baryon to 
entropy ratio rj which in turn can be related simply to Q,bh^ where Jlf, is the ratio of the baryon density to 
the critical density || |6) l^b/i^ _ 3.64 x lO^j^. 

Big Bang nucleosynthesis provides a precise estimate of the baryon density of the Universe by comparing 
to observations of the abundance of deuterium. The most recent observations ]37| constrain the baryon 
density inferred from standard big bang nucleosynthesis to the interval 



0.015 < rJb/i^ < 0.024 (3.27) 

with rjf, being the baryon density expressed as a fraction of the critical density and h has been introduced 
in eqn. ( |2.2| ). In the past year, however, combined results on the CMB anisotropics from the balloon 
experiments BOOMERANG and MAXIMA and observations from type la supernovae have set new limits 
for the baryon density (independent from the observation of abundances) 

0.022 < Vlbh^ < 0.039 (3.28) 

Although these results may be preliminary and require other observations for definite confirmation, taken 
together they signal the possibility of non-standard nucleosynthesis. These new observational constraints 
on fli, have rekindled the interest on the possibility of inhomogeneous nucleosynthesis with encourag- 
ing results for the abundances of deuterium and helium but the abundance predicted by these models of 
inhomogeneous nucleosynthesis for heavier metals is still well below the observed values. 

The important point in this discussion and the main reason for delving into the subject of nucleosynthesis 
is to highlight two important factors that can be infiucnced by the QCD phase transition: i) the helium 
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abundance depends crucially on the neutron to proton ratio, although it is only slightly sensitive to the 
baryon to photon ratio, ii) The abundance of heavier elements up to Li depends in an important manner 
on the baryon to photon ratio. The point being made here is that if the QCD phase transition is of first 
order, there is the potential for producing inhomogeneities that modify the neutron to proton ratio, hence the 
helium abundance, as well as the baryon to photon ratio, and therefore the abundance of heavier elements. 
Hence if the value of ili,h^ favored by BOOMERANG-MAXIMA and Type la supernovae searches stand up 
to further and deeper scrutiny, and are still in discrepancy with standard big bang nucleosynthesis perhaps 
inhomogeneous nucleosynthesis at the QCD phase transition could provide an explanation or be a part of 
the explanation. 

3.5 Recombination, LSS and CMB 

Continuing down the ladder of energy scales, the next step corresponds to energies of order eV determined 
by the binding energy of the hydrogen atom. Actually two different (but related) processes occur at this 
energy scale: a) (re) combination of electrons and protons into neutral hydrogen atoms which is described 
by the chemical reaction 



e~+p^H + 'y (3.29) 

b) photon decoupling (or freeze out): electrons and photons are in LTE through Thompson scattering 
+ 7 ^ e~ + 7. This electromagnetic process has the (large) cross section 



ar - ^ = 6.6 x IQ-^^cm^ (3.30) 

When the reaction rate = u^ut with the electron density becomes smaller than the expansion 
rate (or alternatively the photon mean free path becomes larger than the Hubble radius) the photons do 
not scatter any more, their distribution freezes and redshifts with the expansion. When the free electrons 
combine with protons binding into neutral hydrogen atoms the reaction rate becomes very small. Detailed 
calculations |^ reveal that recombination and decoupling occur at a temperature ~ iQQQK at a redshift 
z ~ 1100. The reason that the recombination temperature is not given by the binding energy of the hydrogen 
atom (= 13.6 eV) is again a consequence of the fact that the baryon to photon ratio is so small as given 



by (3.17) that even a small number of energetic photons in the blackbody tail can photodisintegrate the 
hydrogen atoms formed at T > Tn. 

The important point here is that at a redshift z ^ 1100 and temperature Tfj ~ 3000 K photons freeze-out, 
i.e, they do not undergo any further scattering. This occurs at a time tj^ ^ 300000 years after the Big Bang 
which defines the Last Scattering Surface (LSS). Photons have been travelling to us from the last scattering 
surface without scattering for the last ^15 billion years and carry information of the LSS. 

The discovery of the temperature anisotropics in the CMB by the COBE satellite with AT/T ~ 10"^ on 
angular scales from 90° down to about 2° raises two fundamental questions: 

• The Horizon Problem: photon decoupling occurs at t ^ 300000 years after the Big Bang at a 
redshift of 2: ~ 1100. At this time the Hubble radius is dnitR) ^ 100 Kpc which determines the size 
of a correlated region. Today the size of this correlated patch is ~ (1 -I- z)dH{tR) ~ 100 Mpc, however 
the Hubble radius today is ^ 3000 Mpc. Therefore we would expect that in the Hubble volume today 
there are about ^ 30000 of these regions. Inside each of these regions the temperature fluctuations 
would be correlated but the regions would be completely uncorrelated between them. Each correlated 
patch today subtends an angle of about 2° which is the angular scale of the LSS today. However the 
CMB is homogeneous to 1 part in 10^ on angular scales between 2° up to 90° as revealed by COBE. 
Thus the question: how did these regions that were outside the causal horizon of each other at photon 
decoupling manage to establish correlations and be so homogeneous in temperature?. 
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• The origin of the temperature anisotropies: what "seeds" the temperature inhomogeneities?. 
Simple statistical fluctuations are far too small. This can be seen by considering the simple case of 
a lump of matter with about one solar mass which has about lO^'' particles. Statistical fluctuations 
will be of order ^ 10~^^ which cannot be reconciled with a AT/T ^ 10^^. Obviously for much larger 
masses the fluctuations are much smaller. 

Inflation provides a natural answer to both questions: the temperature anisotropies are caused by small 
inhomogeneities in the matter density which originated in the quantum fluctuations during the inflationary 



era as described in section 3.1.1 above. These density perturbations are given by eqn. (3.7) with physical 
wavelengths that are deep inside the Hubble radius and therefore were in causal contact during inflation 
(see fig.(|l|)). When these wavelengths cross the horizon (becoming super-horizon) during inflation, they are 
no longer affected by microphysical processes, but carry the initial correlations that they had while inside 
the Hubble radius. When they re-enter the Hubble radius right before the LSS, these density fluctuations 
provide anisotropies in the gravitational potential. As the photons fall in these anisotropic potential wells 
their wavelengths are redshifted, regions with larger densities provide larger gravitational potential and 
photons are redshifted more and are therefore cooler. These regions were in causal contact during inflation, 
became causally disconnected from any physical process until they re-entered again the Hubble radius, 
carrying the original correlations. Although there are many details in the relation between temperature and 
density anisotropies |3| this simple argument illustrates some of the most important physics for the 

anisotropies on the scales measured by COBE. A detailed analysis ||, Q || yields 



-Tf o \lSS (•3.31) 

To 3 p 

The main point of this discussion and the reason for delving on this subject, is to highlight that the CMB 
reveals information on processes that occurred very early, for example during inflation, when the physical 
wavelengths of the perturbations that generate the temperature anisotropy first crossed the horizon, or very 
late, near recombination right after they re-entered. To precisely highlight this point we have purposely 
studied in section 3.1 the influence of a phase transition on the power spectrum of density perturbations, 
and therefore by eqn. (3.31) on the temperature fluctuations. That is to say, and we emphasize, that if 
the CMB reveals any observable consequences of phase transitions-the main theme of this lecture- these 
are not part of the standard model: either these phase transitions ocurred at the time of inflation (perhaps 
some GUT?) or at a very low energy scale of 0{eV). The acoustic peaks in the multipole expansion of the 
temperature anisotropies can give indirect information on possible phase transitions through constraints on 
ilb/i^ and the cosmological parameters H. 



3.6 Galaxy formation and on to Stars... 

The temperature anisotropies measured in the CMB give the clue to the origin of large scale structure and the 
formation of galaxies. As described above, the temperature fluctuations of the CMB reflect fluctuations in 
the density of matter. The currently accepted scenario for the formation of galaxies and large scale structure 
begins with the growth of small (linear) perturbations in the density under the relentless action of gravity. 
To understand the main ideas in a simplified manner, consider a small fiuctuation Sp in a background of 
constant density pQ. A region of larger density has a larger gravitational potential which in turn causes 
this region to become denser, i.e, gravity tends to make this overdense region even more overdense. The 
dynamical time scale for the free-fall collapse of this region is 



tff ^ - — (3.32) 

VGnPo H 

with Gn = ^/Mpi being Newton's gravitational constant. However a change in density induces a proportional 
change in pressure which tends to restore hydrostatic equilibrium. The proportionality constant is determined 
by the speed of sound. 
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Sp = c^Jp 



(3.33) 



Consider a perturbation of wavelength A, if the time it takes a pressure wave to restore hydrostatic 
equihbrium over this wavelength is smaller than the free-fall time (3.32), then the restoring force from the 
pressure wave prevents the gravitational collapse. On the other hand if the free-fall time for collapse is 
shorter than the time scale for the pressure wave to restore equilibrium over the distance A, gravitational 
collapse is unhindered. The balance between these two time scales determines the Jeans wavelength 



A ■ [ ^ collapse , . 

^ \/GnPo ' 1 < no collapse 

this is the Jeans instability. In a non-expanding geometry this instability leads to the exponential growth 
of density perturbations for A > Aj, but in an expanding geometry the instability has to catch up with the 
expansion and as a result small perturbations grow as a power law with time . In particular for a matter 
dominated Universe (as is the case after recombination) the density profile for A > Aj grows in time as 



Po 

or equivalently 

Spjt) 
Po 



a{t) cx t^ (3.35) 



(l + z)-' (3.36) 



To get an idea of the order of magnitude of the Jeans scale it is illuminating to compute the Jeans mass 
(mass contained in a volume of radius Aj/2) at the time of recombination. At this time most of the matter 
is in the form of hydrogen with an (adiabatic) speed of sound = bkuT /imH leading to |^ 

= 1.6x10^^^ (3.37) 

6 {^oh^)^ 

which is the size of a typical globular cluster. 

Since at the LSS 5p/ po = 3AT/T ~ 10~^ at a redshift zlss ^ 1100, the argument above suggests that 
density perturbations become non-linear (their amplitude becomes of C(l)) at redshift of C(l), which is 
consistent with observations of star- forming regions |3^. Thus the temperature fluctuations in the CMB 
give direct confirmation of the paradigm of structure formation based on the gravitational instability of 
primordial density perturbations. 

The main points of the analysis provided in this section are: i) to argue that the primordial density 
fluctuations that determine the temperature fluctuations of the CMB through ( 3.3l| ) provide the seeds for 



structure formation. Any information from a primordial phase transition that is imprinted in the anisotropies 
of the CMB is also imprinted in the initial power spectrum of density perturbations that seed large scale 
structure formation, ii) to highlight that when the pressure vanishes, gravitational collapse on large scales 
goes unhindered. This last observation relates to the previous discussion of the formation of primordial black 
holes during the QCD phase transition, which if first order implies an anomalously small (if not vanishing) 
speed of sound. 

3.7 Stellar Evolution 101: QGP at the core of Pulsars? 

The basic Jeans gravitational instability is also at the heart of star formation. Protostars condense from 
giant molecular (hydrogen) gas clouds [Eol . A solar mass protostar condenses in typical star forming regions 
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with temperatures T ^ K and densities > 10~^^ gjcriT' the gravitational contraction of the cloud heats up 
the gas until the temperature reaches ~ 10^ K when thermonuclear fusion of hydrogen into helium generates 
enough energy and pressure to halt the collapse. The threshold mass for a star to undergo thermonuclear 
fusion is about 0.08 M©. Stars with smaller masses do not attain the temperature to ignite hydrogen and 
the gravitational pull is counterbalanced by electron degeneracy, these "failed" stars are called brown dwarfs. 
Solar mass type stars live on the main sequence for about 10 billion years and spend most of their lives burning 
hydrogen into helium mainly through the three branches of the pp (proton-proton) chain and more massive 
stars also produce helium through the C(arbon)N(itrogen)0(xygen) chains |Q. When hydrogen is depleted, 
gravitational contraction begins again until the temperature reaches about lO^K when the thermonuclear 
fusion of helium begins at the core while a shell of hydrogen still burns outside the core. This helium 
burning stage produces carbon (through the miracle of the triple alpha process) and oxygen. Computer 
simulations |^0| show that stars with masses up to ~ 8Mq burn hydrogen and helium but the temperatures 
attained at the core are not high enough to ignite carbon. These stars (the sun is one of them) end up their 
life cycles as white dwarfs, with a degenerate CO core, the gravitational pull is counterbalanced by the Fermi 
degeneracy pressure of electrons. Typical radii and densities of these stars are R ~ 10"* Km; p ~ 10^ gr/crrv^. 
The maximum mass for a stellar object to be in hydrostatic equilibrium balanced by Fermi degeneracy 
pressure is given by the Chandrasekhar limit ~ 1.4 Mq. When the density is such that the Chandrasekhar 
limit is reached, electrons at the top of their Fermi surface become ultrarelativistic and the equation of 
state is softer. Stars with masses > 8 — 10 Mq evolve through all the stages of nuclear burning beginning 
with hydrogen, continuing with helium, carbon, neon, oxygen and silicon burning (photodisintegration at 
T ^ 3.5 X 10^ K) which ends in the iron group elements with an ^^Fe core and concentric shells ("onion 
structure") of silicon, neon, oxygen, carbon, helium and hydrogen. Since Fe has the largest binding energy 
per baryon, thermonuclear fusion cannot proceed further. For stars with masses > IOMq the iron core 
reaches the Chandrasekhar limit for support by electron degeneracy and is on the limit of gravitational 
collapse. There are two main factors that trigger the collapse |^0[ i) the temperature of the core is 

now ~ 8 X lO^K enough to photodisintegrate iron group elements by the reaction 



7 + ^''Fe <^ WHe + in (3.38) 

ii) the Fermi energy of the degenerate electrons becomes larger than the neutron-proton mass difference 
1.3Mev and is therefore large enough for electron capture by protons 



e + p ^ n + (3.39) 

this is the process of neutronization. As the electrons that support the gravitational pull with their Fermi 
pressure are captured, hydrostatic equilibrium falters and the core begins to collapse on a dynamical (free fall) 
time scale given by eqn. ( 3.32| ). With the density of the core ~ 10^ g/cnr' the free fall time scale for collapse 



is ~ 1msec. As the collapsed core reaches nuclear matter density, it becomes highly incompressible and its 
equation of state stiffens. The infalling matter from the outside shells bounces from the incompressible core 
and a shock wave is formed. This shock blows off the outer layers resulting in a Supernova (core-collapse 
or Type II) explosion. The details of the explosion and the shock are very complex and involve neutrino 
transport, convection etc. For more details the reader is referred to [^0[ 

What happens after this depends on the mass of the progenitors. Current understanding based on 
numerical evolution |4^, ^ Q suggests that for stars with masses up to about 3GM0 the final result is a 
neutron star, while for masses larger than this the collapse probably leads to the formation of a black hole. 
When the supernova explosion leads to a neutron star, an important feature with observational consequences 
is that during the process of neutronization of the iron core with the Chandrasekhar mass ~ 1.4Mq most of 
the 10^^ protons present in the core are converted to neutrons and ^ 10^^ neutrinos each of about lOMev 
are released and carry most (up to 99%) of the energy ^ 10^'^ergs. When the density of the collapsing core 
reaches about 4 x 10^^ g/cm^, neutrons begin to drip out of the nuclei, this is the neutron drip line and the 
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collapsing core becomes a gigantic neutron. Typical neutron stars have masses in the range 1.4 — 2M0 with 
typical radii ~ 10 km. The most direct evidence for neutron stars are pulsars which are believed to be highly 
magnetized {B ^ IQ^'^G) rotating neutron stars Pulsar periods range from about a millisecond up 

to a second. 

The equation of state for neutron stars for densities below nuclear matter p„ ~ 2 x IQ^^g/cnv' is fairly 
well understood p5| |4^ , while for densities above p„ there is the possibility of hyperon rich matter, pion 
and or kaon condensates, muons, and other exotica ^ ^ Modern theories of superdense nuclear 
matter 43, 4^, predict a composition for neutron stars that is depicted qualitatively in fig.(^ below. 




nuclei and electrons 
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Figure 2: Neutron star composition. 



A crust of i^^Fe) nuclei, electrons and neutrons, followed by an internal region of neutron and proton 
superfluids with the possibility of pion and kaon condensates and hyperons. But more importantly for this 
discussion, current ideas suggest the possibility of a mixed phase of hadrons and quarks, characterized by 
geometric structures (a consequence of the competition between nuclear and Coulomb forces) and ultimately 
a core of deconfined quarks and gluons with T < 1 Mev and p ^ 3 — 5p„, i.e, the core is conjectured to be a 
cold, dense Quark Gluon Plasma. 

Recently the interesting possibility of novel color superconducting phases of cold and dense QCD has 
been proposed and we refer the reader to several review articles on the subject for details on these fascinating 



aspects |47 



In a very tantalizing recent article Glendcnning and Weber proposed that a phase transition to 
quark matter can provide a potential explanation for the anomalous distribution of pulsar frequencies in 
low mass X-ray binaries (LMXBs). Observations of the pulsar frequencies in LMXBs by the Rossi X-ray 
Timing Explorer reveal a spike at a frequency ~ 300 Hz. We will describe this fascinating possibility below 
along with the current experimental efforts to map the QCD phase diagram with ultrarelativistic heavy ion 
collisions. 

The physics is different for Type la Supernovae that are currently at the center of the discussion of their 
use as standard candles to study the cosmological constant These are conjectured to be CO white 

dwarfs in a binary system. The white dwarf accretes mass from its companion until its own mass becomes 
just below (a few percent) of the Chandresekhar limit. Current understanding suggests that at this point 
the CO core begins thermonuclear burning under degeneracy conditions. This results in runaway burning 
since in the Fermi degenerate situation the pressure is almost independent of the temperature (just as in 
the helium flash). The temperature increase at almost constant pressure leads to a runaway thermonuclear 
burning that results in the explosion of the star |^ . What makes these supernovae so special is that their 
composition and mass at the time of the explosion is the same for all Type la Supernovae and this feature 
may explain their similar light curves KM . 
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3.8 Executive Summary: observational consequences of cosmological phase tran- 
sitions? 

The focus of the review of big bang cosmology as well as the astrophysics of (compact) stars is to highlight 
where and when there could be observational consequences of the phase transitions predicted by the standard 
model of particle physics. The picture that emerges from this brief tour through the early Universe, is 
that while current CMB measurements and large scale surveys have the potential for revealing evidence 
for inflationary phase transitions the observable evidence for a standard model phase transition is at best 
indirect. I have argued that inflationary phase transitions, perhaps taking place at the GUT scale could 
leave an imprint in the power spectrum of temperature anisotropics, for example a red tilt as a consequence 
of spinodal instabilities |^] . Or as argued in ref. a step-like feature in the power spectrum of temperature 
anisotropics or in the power spectrum of galaxy clusters may be produced by an inflationary phase transition 
in a supersymmctric theory. Either proposal however, lies beyond the standard model and at energy scales 
unlikely to be probed by any current or future accelerator. The main "problem" with detecting phase 
transitions with the CMB or large scale surveys is that the fluctuations that can seed temperature and 
density inhomogeneities were generated very early, during inflation, i.e, at an energy scale ~ 10^^ Gev, or 
very late during recombination, i.e, at an energy scale ^ 0.3 eV. As is clear from fig.(|^) scales of cosmological 
relevance today, say from ^ lOMpc up to the Hubble radius ~ SOOOMpc, were super-horizon, hence decoupled 
from the microphysics at the time of the EW or QCD phase transitions. 

The indirect evidence for standard model phase transitions is contained in the baryon to photon ratio, 
i.e, the baryon asymmetry and in the possibility of primordial black holes or inhomogeneous nucleosynthesis 
caused by the QCD phase transition. As argued above, it now seems clear that the standard Electroweak 
theory cannot describe consistently baryogenesis because of the smallness of the CP violating parameter(s) 
as well as the large Higgs mass. This analysis leaves the responsibility for observational consequences of 
standard model phase transitions squarely on the shoulders of QCD!. A strong first order phase transition 
from a QGP to a hadron phase can lead to the formation of primordial black holes because the speed of 
sound becomes anomalously small thus gravitational collapse of solar mass clumps is unhindered. There is 
the possibility of inhomogeneities in the baryon density leading to inhomogeneous nucleosynthesis with the 
tantalizing possibility of explaining the new bounds on il^h^ from Boomerang and Maxima. There is also 
the possibility for the formation of strangelets or strange quark nuggets that could be a component of the 
cold dark matter. As described above most of these possibilities rely on the details of the QCD FT and the 
parameters of QCD that determine the main dynamical aspects of the transition. While the lattice gauge 
theories program has the potential for computing many of these quantities, it is clearly important to study 
experimentally the phase transitions of the standard model. As mentioned above the direct study of the 
EWPT is certainly out of the reach of current and most certainly near future accelerators: to achieve the PT 
temperature T ^ 100 Gev requires an energy density e ^ lO^^pn- The QCD phase transition is perhaps the 
only phase transition of the standard model that can and is being studied with accelerators. At the QCD PT 
temperature T ~ 200 Mev the energy density required is e ~ 2 — 4p„ which is the energy density achieved 
in ultrarelativistic heavy ion collisions. I now describe the current (and future) program to study the QCD 
phase transition(s) with ultrarelativistic heavy ion collisions, as well the recent fascinating suggestion that 
there could be observable hints of quark matter in the core of pulsars. 

4 Relativistic heavy ion collisions and pulsars open a window to 
the early Universe: 

The program of relativistic heavy ion collisions whose primary goal is to study the phase diagram of QCD 
began almost two decades ago with the fixed target heavy ion programs at the AGS at Brookhaven and 
the SPS at CERN. Recently |3^ a summary of the results of these efforts mainly through the Pb + Pb 
experiments at SPS-CERN provided very exciting "evidence" in favor of the existence of the QGP. As 
impressive as this body of evidence is, the consensus in the field is that while the evidence is very suggestive 
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it is far from conclusive. As SPS shuts down to pave the way for the forthcoming Large Hadron CoUider 
wherein the Ahce program wiU continue the search for the QGP, the torch now passes to the R(elativistic) 
H(eavy)I(on)C(omder) at BNL. 



4.1 RHIC and LHC seek the QGP: the big picture 

The Relativistic Heavy Ion Colhder (RHIC) at BNL is currently studying Au + Au collisions with center of 
mass energy ^/s ^ 200AGev and luminosity ~ lQ^^cm~^s~^ . The future ALICE heavy ion program at LHC 
is expected to study Pb+Pb collisions with cm. energies up to ^/s ^ 5ATev and luminosities ~ lO^^cm^^s"^. 
In these collisions the heavy nuclei can be pictured in the CoM frame as two Lorentz contracted pancakes, 
for example for Au + Au collisions the size of each "pancake" in the direction transverse to the beam axis 
is about 7 fm. At RHIC and LHC energies most of the baryons are expected to be carried away by the 
receding pancakes (the fragmentation region) while in the region of the collision a large energy (density) is 
deposited in the form of quark pairs and gluons. At least two important mechanisms for energy deposition 
in the collision region are at work i) the establishment of a strong color electric field (flux 

tube) that eventually breaks up into quark-antiquark pairs when the energy in the field is larger than the 
pair production threshold and ii) the partons (quarks and gluons) inside the colliding nuclei interact and 
redistribute their energy in a "parton cascade" iQ. An estimate of the energy deposited in the collision 
region has been provided by Bjorken p5| 



1 dEr 

^ = — ^~r~ (4-1 

TQirRj^ dy 

with To ^ Ifm/c, Ra ^ 7fm for Au and dET/dy is the transverse energy per unit rapidity which is measured. 
At RHIC the transverse energy per unit rapidity is expected to be in the range 500 — 900Gev, giving 
for the energy density e ^ 4 — 6Gev/fm^ translating this value into the temperature of a quark-gluon 
gas leads to T > 400Mev larger than the expected critical temperature Tc ^ 160Mev. The evolution of 
the parton distribution functions reveals that thermalization of quarks and gluons occurs on time scales 
0(0.3 — 0.7)fm/c with gluons thermalizing first jsj, After the quark-gluon plasma achieved LTE, the 
evolution is conjectured to be described by hydrodynamic expansion [^6| . As the QGP expands and cools, 
the temperature falls near the critical temperature and the confinement and chiral phase transitions occur. 
Upon further cooling the quark-gluon plasma hadronizes, the hadrons rescatter until the hadron gas is dilute 
enough that the mean free path is larger than the mean distance between hadrons. At this point hadrons 
"freeze-out" and stream out freely to the detectors from this last scattering or freeze-out surface. This picture 
is summarized in fig. (^ below. 

Current estimates based on this picture and on detailed numerical evolution suggest that at RHIC the 
quark gluon plasma lifetime is of order ~ 10/m/c while the total evolution until freeze-out is ^ 50— 100/m/c. 



4.1.1 Hydro, LGT and the EoS: 

Although the evolution of the quark gluon plasma from the initial state described by the parton distribution of 
the colliding nuclei until freeze-out of hadrons clearly requires a non-equilibrium description, a hydrodynamic 
picture of the evolution is both useful and experimentally relevant |^6|. In the approximation in which 
quarks and gluons are strongly coupled in the sense that their mean free paths are much smaller than the 
typical wavelength for the variation of collective phenomena, the QGP can be described as a fluid in LTE. 
Hydrodynamics results from the conservation laws applied to the fluid form of the energy momentum tensor 
and the conserved currents. In particular for a fluid in LTE, the energy momentum tensor is of the form 

T''" = {e+p)u''u'' -pg'"' (4.2) 

with e,p the energy density and pressure respectively, g'^'^ is the metric and 
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Figure 3: Heavy Ion Collisions. 



= 7(x)(l,i/(x)) (4.3) 

is a local 4- velocity vector with j{x) the local Lorentz contraction factor. In this description the dynamical 
evolution of the fluid is obtained from the conservation laws of energy-momentum, baryon number and 
entropy ^ and to close the set of equations one needs an equation of state p — p{e). 

In order to make a definite connection with the cosmological setting as well as for it simplicity and 
usefulness to describe qualitatively the expansion and cooling of the QGP we now briefly summarize Bjorken's 
model for longitudinal expansion [^3| |55| , p6| . There are three important ingredients in Bjorken's models: 
i) the "central rapidity" region where the energy deposited and the formation of the QGP takes place is 
well separated from the "fragmentation region" which is the region of the receding pancakes where most of 
the baryons are (see fig.(^). ii) The thermodynamic variables are invariant under boosts along the beam 
(longitudinal) axis, this is based on the observation that the particle distributions are invariant under these 
boosts in the "central rapidity region" and iii) the expansion only occurs along the beam axis, i.e, longitudinal 
expansion. In the baryon free region, local thermodynamics implies that the entropy density in the local 
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rest frame of the fluid is related to the energy density and pressure as 

e + p = Ts 

It is convenient to introduce the space time rapidity variable rj and proper-time r 



(4.4) 



^=-ln 



t + z 
t - z 



(4.5) 



with z the coordinate along the beam axis (longitudinal). Introducing the fluid rapidity Q by writing the 
local velocity (4.3) as 



= (cosh(0), 0, 0, sinh(0)) (4.6) 

Furthermore Bjorken's model assumes the fluid to be composed of free streaming particles for which 
Vz — z/t in which case the fluid rapidity 9 becomes the space-time rapidity 77. Boost invariance along the 
longitudinal direction entails that e,p,s,T are all functions of proper time only (here s,T) are the entropy 
density and temperature ||5^, |5^ . 

The energy-momentum conservation equation d^T^'^ leads to two equations by projecting along the 
direction and perpendicular to it using the projector g'"' — u^^u^ . Under the assumption that the central 
region is baryon free, the conservation of entropy and energy and momentum lead to the following equations 
(for details see g |l], ||, ||) 



de{T) 

dT 
ds{T) 

dT 



-(e+p) = 

T 

r 



(4.7) 
(4.8) 



Assuming an equation of state p{t) = p{e{T)) and combining eqns. (4.7 4.81 ) with the thermodynamic 
relation (4.4) (for baryon free plasmas) one finds that the evolution of the temperature is given by 



T dT 



dP 

de 



(4.9) 



Which for constant speed of sound results in the cooling law 



r(T) = ro(^)'^ (4.10) 



T 



Obviously the form of the equation (4.7) is similar to the energy conservation equation in a homogeneous 
and isotropic metric given by (2.8) with the expansion rate a/a — l/3r when the proper time r is identified 
with the comoving time t in the cosmological setting. The similarity becomes even more remarkable for the 
case of the QGP being modelled as a radiation fluid (which is expected to be a good approximation at high 
temperature) since in this case ~ 1/3 and the connection with a radiation dominated cosmology with 
scale factor a{t) (x ti is evident. 

To find the general evolution equations for the QGP under the conditions specified above, an equation 
of state is needed. 

It is at this stage where the connection with the lattice gauge theory (LGT) program is made. LGT 
obtains the thermodynamic functions in equilibrium and these are input in the hydrodynamic description as 



local functions of space-time under the assumption of LTE. Fig. (|4|) below is an example of recent results 
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Figure 4: e/T^ and p/T"*, vs T/Tc for a 16^ x 4 lattice.From |p 



from LGT for the energy density and pressure (both divided by to compare to a free gas of massless 
quarks and or gluons) as a function of T/Tc with Tc = 170 Mev. 

It is clear from these figures that there is a sharp decrease in the energy density and pressure ai T — 
the value of the energy density at Tc is, from the figure Sc ^ 6T^. Furthermore the large temperature 
behavior is not quite given by the Stephan-Boltzmann law ((SB) in the figures), suggesting that even at large 
temperatures the plasma is not described by free quarks and gluons up to temperatures T ~ 4,.0Tc 700 Mev. 
This discrepancy is currently explained in terms of quasiparticles rather than free particles (see below). 
Nevertheless, what transpires clearly from these figures is that for temperatures T > Tc the plasma is 
approximately a gas of almost free massless quarks and gluons, while for T «Tc the plasma is dominated 
by heavy particles. This is the picture that deconfinement-confinement phase transition suggests. 

A hydrodynamic description combined with the non-perturbative lattice gauge theory program to extract 
the EoS of the QGP provide a basic quantitatively reliable picture of the dynamical evolution of the QGP. 
The approximation of ideal (or perfect) hydrodynamics can be relaxed by considering viscous hydrodynamics 
by including the bulk and shear viscosity in the energy momentum tensor |59|| . 

4.2 Thermalization, quasiparticles and the EoS: 

The initial stages of evolution of the QGP beginning from the distribution of partons (quarks and gluons) 
in the colliding nuclei is studied with a semiclassical transport approach that includes perturbative QCD 
cross sections with screening corrections [^^. These studies reveal that gluons and quarks thermalize on 
time scales < 1/m/c with gluons thermalizing first. While a classical (or semiclassical) transport approach 
may not be completely justified in this regime, a qualitative and quantitative idea of the relaxation time 
scales can be obtained by assuming a thermalized state and studying how a small departure from equilibrium 
relaxes. The relaxation rate is given by F = ncr with n the particle density and a a typical scattering cross 
section. Consider the relaxation of quarks in a bath of gluons at temperature T, one gluon exchange yields 
a typical cross section 



y2 



(4.11) 



with tts = 5^/471 and g is the gauge coupling constant and we have assumed that the typical energy of the 
exchanged gluon is ^ T. For T >> m with m the typical masses, n ocT^ leading to 



a 



(4.12) 
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Gluons of typical energy T are hard and lead to smaller cross sections, while the exchange of soft gluons 
leads to larger cross sections and is more efficient for thermalization. However soft gluons are screened by 
the medium. The estimate above, eqn. (4.11) is reliable for hard gluon exchange, but overlooks the fact 
that soft longitudinal gluons are Debye screened. A careful calculation reveals that the Debye screening 

"^T^. Furthermore, 
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length mz) gT therefore the denominator in eqn. (4.11) must be replaced by T 
while longitudinal (electric) gluons are Debye screened, transverse (magnetic) gluons are only dynamically 
screened by Landau damping |6l|, Including these screening effects it is found that the quark relaxation 
rate is given by 



r~a,rin(^-j (4.13) 

for Us ^ 0.3; T ~ SOOMev one finds the typical relaxation time T^ei — ^ 0.8 — Ifm/c. 

The main point in bringing up this estimate is to highlight the fact that in a medium there are many-body 
effects that screen or "dress" the particles, and as a result, even for high temperatures the plasma is not 
described by a free gas of quarks and gluons, but in terms of "dressed" weakly interacting quasiparticles. In 
a thermalized QGP these quasiparticles have a typical "thermal mass" of order gT |Q. A framework to 
include these screening corrections consistently in the perturbative expansion has been put forth in ref. 
and is referred to as the hard thermal loop (HTL) program. In fact the energy density and pressure obtained 
by the LGT program depicted in fig. (jj) have been recently reproduced analytically in terms of thermal 
quasiparticles with a typical mass of order gT |l63| . Thus the combination of numerical results from LGT and 
the analytic resummation provided by the HTL program lead to a picture of the QGP at high temperatures 
in terms of weakly interacting quasiparticles with thermal masses gT. 



4.3 Predictions and observations: 

Several experimental signatures had been associated with the formation of the QGP and the heavy 

ion programs had focused on several of them. Here I summarize the observables that have been proposed 
as telltales of a QGP and the data gathered by the SPS from Pb + Pb and Pb + Au collisions, which taken 
together provide a hint of evidence ^ ^ for the effects associated with a QGP (although many if not 
all of them could have alternative explanations). 



4.3.1 J/^f suppression: 

The J/^* is a cc bound state that is very narrow and that if produced in the early stages of the collision 
can probe the QGP because its lifetime is longer than that of the QGP and decays into dilepton pairs which 
leave the plasma without scattering. The original suggestion is that when the screening length of the 
color force is smaller than the size of the cc bound state, this narrow resonance will "melt" . This original 
argument thus suggests that a suppression of the J/^* yield could provide information on the existence of 
the QGP. Alternatively a suppression could result from scattering with hard gluons in the plasma and the 
"dissociation" of the bound state |]65|| . A normal suppression of charmonium is expected on the grounds that 
once formed the cc bound state interact with other nucleons inside the nucleus. This expected suppression 
is studied in proton nucleon collisions and extrapolated to nucleon-nucleon collisions. This is considered the 
"normal" suppression in contrast to the "abnormal" suppression expected from the presence of a plasma. 
Fig. (|) below shows the data gathered by the NA50 coUaboration at the SPS-CERN ||6|. 

This striking figure clearly reveals an "abnormal" suppression when the energy density is e > 2.5Gev/fm'^. 
The energy density in this figure has been computed with Bjorken's formula ( [4.l| ). Further analysis jsj, |3^ 
reveals that this anomalous suppression cannot be accounted for by coUisional dissociation through final 
state hadronic interactions. 

A recent report from the NA50 collaboration at CERN-SPS ||6^ presents the combined data for J/'^ 
suppression from the NA38 and NA50 experiments. The analysis of the data reveals that while for the most 
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Figure 5: Anomalous J/?/; suppression as a function of the initial energy density. From 




peripheral (largest impact parameter) collisions the suppression can be accounted for by nuclear absorption, 
there is no saturation in the suppression in the most central Pb + Pb collisions and that the observed 
suppression pattern can be naturally understood in a deconfinement scenario. This report concludes that 
the J/'^ suppression pattern observed in the NA50 data provides significant evidence for deconfinement of 
quarks and gluons in Pb + Pb collisions. 

4.3.2 Electromagnetic probes: dileptons and direct photons 

Electromagnetic probes: e+e^ or fJ.'^ fJ-^ dilepton pairs and direct (prompt) photons are prime probes of the 
hot plasma | |68| , since once produced they leave the plasma without further interactions because their mean 
free path is much larger than the typical size of the plasma. 
Dileptons: 

The p vector meson has particular relevance in this regard because it decays into dileptons and its lifetime 
is ~ 1/m/c, therefore once it is produced in the hadron gas it decays within the hot hadronic plasma and 
the produced dileptons carry directly this information. Thus while dileptons produced from the decay of the 
p meson do not yield evidence of the earlier stages in the QGP, they do nevertheless offer information on 
the hadronic stage. Figure (^ below presents the data gathered by the CERES-NA45 collaboration for the 
invariant mass spectrum for electron-positron pairs from 158 AGev Pb + Au collisions at the SPS-CERN. 

The solid line represents the expected spectrum from the decays of hadrons produced in proton-nucleon 
and proton-proton collisions extrapolated to Pb + Au collisions and is the sum of the contributions shown in 
the graph. There are two remarkable features in this graph: a clear enhancement of dileptons in the region 
250 Mev < Mg+e- < 700 Mev and that instead of the p meson peak at nip = 770 Mev there is a broad 



28 




Figure 6: Invariant mass spectrum of e+e pairs from Pb + Au collisions at 158 AGev at the SPS-CERN 
(CERES/NA45). From |6|. 



distribution. The excess of dileptons in the small invariant mass region cannot be explained by charged pion 



annihilation , 35 . What is remarkable in this data is that the dilepton enhancement is below the putative 
p peak and that there is no hint of the p at 770 Mev!. The current understanding of these features is that 
the medium effects result in a shift in the p meson mass as well as a change in its width |7^, ^ |5^. Thus 
while this interpretation does not directly yield information on the QGP, it does support the picture of a 
hot gas of hadrons, mainly pious which is the main interaction channel of the p vector meson. 
Direct photons: 

Direct photons are conceptually a clean direct probe of the early stages of the QGP. Photons are produced 
in the QGP by several processes: gluon-to-photon Compton scattering off (anti)quark q{q)g ~* q{q)^ and 
quark-antiquark annihilation to photon and giuon qq — + 57 and to the same order (see Kapusta et. al. 
and Aurenche et. al. in [68| ) (anti)quark bremsstrahlung qq{g) — > qq{g)^ and quark-antiquark annihilation 
with scattering qqq{g) q{g)l- Detailed calculations including screening corrections p8) reveal that direct 
photons from the QGP could provide a signal that could be discriminated against that from the hadronic 
background. This work indicates the theoretical feasibility of direct photons as direct probes of the early 
stages of the QGP. 

Recently the WA98 collaboration at SPS-CERN reported their analysis for the first observation of direct 
photons from Pb + Pb collisions with y/s — IbSAGev |^^. Their data is summarized in the fig.(0) below. 

The transverse momentum distribution of direct photons is determined on a statistical basis and compared 
to the background photon yield predicted from a calculation of the radiative decays of hadrons. The most 
interesting result is that a significant excess of direct photons beyond that expected from proton-induced 
reaction at the same -^/s is observed in the range of transverse momentum greater than about 1.5 GeV/c in 
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Figure 7: Ratio of total measured yield of direct photons to hadronic background vs. transverse momentum 
from Pb + Pb at 158 AGev for peripheral and central collisions. From [[72|. 



central collisions. 

A detailed analysis of the data and the theoretical expectations based on the calculations of direct photons 
from an equilibrated QGP was recently performed |7^. The conclusions of that analysis is that while it is 
not clear if SPS has reached the energy density to form the QGP, the data supports indications of a hot and 
dense phase that could be the precursor of the QGP. 

More recently [73| it has been suggested that non-equilibrium effects in an expanding QGP formed in 
RHIC collisions could lead to an enhancement of direct photons in the region of transverse momentum 

^ 2 Gev. 

4.3.3 Strangeness enhancement 

Strangeness enhancement along with chemical equilibration are some of the earliest proposals for clear sig- 
natures of the formation of a QGP |7^. The main idea is based on the estimate that the strangeness 
equilibration time in a hot QGP is of the same order as the expected lifetime of the QGP 10/m/c) pro- 
duced in nucleus-nucleus collisions. Two important aspects of this estimate make strangeness enhancement 
a prime candidate: if strangeness attains chemical equilibrium in the QGP, this equilibrium value is signif- 
icantly higher than the strangeness production in nuclcon-nucleon collisions. Also strangeness production 
through hadronic rescattering or final state interactions was estimated to be negligibly small |Q. In the 
QGP, color deconfinement leads to a large gluon density that leads to the creation of ss pairs, furthermore 
chiral symmetry makes the strange quark lighter thus lowering the production threshold. This situation 
is in contrast to the case of hadronic rescattering or final state interactions where the production of pairs 



of strange quarks has large thresholds and small cross sections |34, The usual measure of strangeness 
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enhancement is through the ratio 



A, = 



2(ss) 
{uu + dd) 



(4.14) 



Fig. (||) below displays As as a function of ^/s the energy of the collision for nucleon-nucleon as well as 
nucleus-nucleus collisions {S + S, S + Ag, Pb + Pb) at SPS. 
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Figure 8: The ratio A^ as a function of -^5, the energy of the collision for nucleon-nucleon and nucleus-nucleus 
collisions. From IITm 



It is clear that nucleus- nucleus collision is creating a different environment that enhances the formation 
of strangeness. As mentioned above, and confirmed by detailed numerical evolution of kinetic models p6|, 
the enhancement cannot be explained by hadronic scattering. Thus the enhancement displayed by fig. ^) 
provides a very exciting hint that nucleus-nucleus collisions at SPS are creating a hot state of matter. 

4.3.4 Collective flow and the Eos: 

After the collision and the formation of the QGP, the pressure in the plasma drives its expansion and cooling. 
Pressure gradients give rise to changes in the velocity distribution and the collective motion of the plasma 
which in turn characterize the dynamical spatial and momentum correlations of the particles in the plasma. 
The collective motion of the plasma driven by the internal pressure is referred to as collective flow. This 
flow provides information on the equation of state (Eos) of the plasma as is manifest in the hydrodynamic 
description described in section 4.1.1 above. 
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The collective radial expansion of the plasma is typically assessed by looking for deviations of momentum 
distributions from thermal although perhaps a more clear assessment is obtained by comparing spectra of 
particles with different masses [ [77| . 

If the transition from the QGP to the hadron gas phase is first order and occurs in LTE, then the 
two phases coexist at the same temperature and pressure is constant through the transition (through the 
Maxwell construction) . However the ratio of pressure to energy density decreases and reaches a minimum 
at a particular energy density Sgp ~ lAGev / fm^ known as the softest point of the Eos jTSf . As the energy 
density passes through the softest point the effective speed of sound becomes anomalously small, and the 
small pressure cannot accelerate effectively the matter and the flow stalls. There is a host of numerical 
simulations based on dynamical transport approaches that predict noticeable variations in collective flow for 
different equations of state (see for example [78|). 

The measurement of collective flow and an assessment of the EoS is an integral part of the experimental 
program at SPS and RHIC. There are already interesting results that prove the viability of the study of flow 
to understand the dynamical evolution of the hadronic component ||79[| . But while there is an important and 
exciting body of results from SPS and some recent results from STAR at RHIC it is still rather difficult 
and perhaps premature to extract clear information on the nuclear EoS. 

The main point to bear, however, is that the EoS is experimentally accessible through measurements of 
flow. Hence the experimental program in relativistic heavy ion collisions has access to the EoS which as 
described above is a very important aspect of the QCD phase transition. 



4.3.5 Other predictions... 

There are a variety of other "predictions" that purport to describe potentially observable signatures of a 
QGP and we mention here a few that could also be of potential relevance for early Universe cosmology. 

As mentioned above, there are two phase transitions ocurring at about the same temperature, the 
confinement-deconfinement and the chiral symmetry breaking phase transition. Most of the observables 
described above refer to the confinement-deconfinement phase transition. The chiral phase transition refers 
to the breaking of the (approximate) symmetry corresponding to independent rotations of the right and left 
handed components of the u,d quarks. The three pions are the (quasi) Goldstone bosons resulting from the 
breakdown of chiral symmetry. 

If the chiral phase transition occurs out of equilibrium, it is conceivable that strong fluctuations of the pion 
field could emerge. These fluctuations had been given the generic name of "disoriented chiral condensates" (or 
DCC's) and while there could be different manifestations of these DCC's, all bear in common the notion 
of large amplitude, coherent pion "domains" . These large fluctuations would result in large regions in which 
isospin and probably charge are correlated. Furthermore if these domains are produced via long-wavelength 
(spinodal) instabilities, such as those mentioned within the context of a supercooled phase transition in 



section 3.1 above (see the discussion below eqn. (3.12)) there could be a host of observables associated with 



these domains: anomalous distribution of pions at small transverse momentum, enhanced production of 
direct photons with a distinct polarization asymmetry (net helicity) and strong fluctuations in the particle 
multiplicity on an event by event basis | |8l[ | . All of these are potentially important phenomena if they occur 
during or after the QCD phase transition in the early Universe. Because the pions are pseudoscalar particles, 
if they are produced coherently through spinodal decomposition and they eventually produce photons, this 
could give rise to primordial magnetic flelds with net helicity because of the polarization asymmetry of the 
produced photons. This possibility, however, must be analyzed further for a reliable assessment. 

There has been a substantial experimental effort by the WA98 collaboration at SPS-CERN and by the 
Minimax collaboration at the Tevatron (Fermilab) to detect the signatures associated with DCC's, but so far 
the search has yielded negative results Furthermore, recently the NA44 collaboration at SPS-CERN has 
reported its experimental study of critical fluctuations on an event-by-event basis in Pb + Pb collisions. The 
analysis of the data does not reveal any large fluctuations that could be associated with critical phenomena. 

These results disappointing as they may be should not be interpreted, yet, as the implausibility of the 
physical mechanism. It is conceivable that the temperature region probed by SPS is not high enough to 
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lead to a strongly out of equilibrium chiral transition, which would be needed for the above mechanism to 
be viable Q. 

Thus the study of critical fluctuations and the possibility of large pion domains will be continued at 
RHIC. 

4.4 Little bang vs. Big Bang: 

Having reviewed the standard Big Bang cosmology and the experimental effort to study the QCD phase 
transition(s) and to map the QCD phase diagram with accelerators, we now compare the settings for the 
QCD phase transition during the Big Bang to that prevailing in accelerator experiments, i.e, the "little 
bang". The QCD phase diagram and the different regions studied with accelerator experiments as well as 
the region of temperatures and chemical potential prevailing during the first /zsec after the Big Bang is 
depicted in figure (for details on the CS (color superconducting) phase (s) see [0). For a thorough 
review of the experimental aspects to study the phase diagram see Pq] . 
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early universe 
I 




[GeV] 

Figure 9: QGP phase diagram. Adapted from |Q 

There are important similarities as well as differences between the situations in which the QCD phase 
transition occurs in the early Universe and in accelerators which are worth summarizing. 

• Space-time scales: The space-time scales are very different: in the early Universe the QCD phase 
transition occured at i ~ 1 — 10 x 10~^ sees after the Big Bang when the size of the horizon was 
dn '^-^ 10 Km the expansion rate was H 10^ secs~^. Quark and gluon mean free paths are of order 
~ 1 fm, hence conditions for LTE prevailed. 

In the present Universe, i.e, in heavy ion collision experiments, the lifetime of the QGP is estimated 
to be 10 fm/c and the typical size is ~ 10/m, the expansion time scale is ^ Ic/fm while the 
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mean free paths are still of order ^ Therefore departures from LTE and non-equilibrium effects 

could play an important role in heavy ion collisions, while in the early Universe LTE is a very good 
approximation. 

Baryon density: In the early Universe, the QCD phase transition ocurred in an almost baryon free 
environment. This is because the entropy was dominated by photons (and neutrinos) with a baryon 
to entropy ratio ns/s ^ 10^^. At RHIC and LHC energies it is expected that the central collision 
(central rapidity) region will be almost baryon free, with most of the baryons in the fragmentation 
regions and the entropy dominated by (almost massless) pions with a very small ratio nB/n-,^ in the 
central region. 

LSS,CMB... freeze out,HBT There is a direct analogy between photon decoupling, the last scat- 
tering surface (LSS) in cosmology and the freeze-out of hadrons: in the case of photon decoupling, 
the mean free path for Thompson scattering becomes of the same order as the Hubble radius and 
the photons consequently free-stream. A similar situation occurs for hadrons, when their mean free 
path becomes larger than the typical separation of hadrons, there is no further rescattering and the 
hadrons free stream towards the detector. Just as in the case of the CMB where correlations in the 
temperature anisotropics averaged over the sky give information on the last scattering surface, there is 
a similar technique for hadrons, mainly pions, which makes use of the Hanbury-Brown-Twiss (HBT) 
interferometric effect 85, Consider a source that emits identical particles from positions Pi 



and P2 and these particles are later observed at points P3 and P4 as envisaged in fig. (|l^). Because 
of the symmetrization of the pion wavefunctions, both emission points contribute to the observable at 
both reception points, even if the particles are non- interacting. The correlations in the momenta of 
the two pions are studied by defining 

C{p,?)- Jft^,^ (4.15) 

with P2(p, 9) the joint probability of two pions with momenta p and q respectively and the Vi are the 
individual probabilities, so that if C(p, q) = 1 the events are uncorrelated. 

The Fourier transform of V2 gives information on the space-time structure of the source that emitted 
the pions [^4| [2l|| . Thus by studying pion interferometry in heavy ion collisions one learns about 
the space-time structure of the hadronic gas at the freeze-out surface, just as the multipole expansion 
of the correlation function of the temperature anisotropics gives information about the cosmological 
parameters at the last scattering surface Q. In particular HBT interferometry can be used to study 
signals from a first order phase transition as well as coherent pion production (DCC's) and is 
therefore an important diagnostic tool for observable consequences of the QCD phase transition(s), 
just as the analysis of the temperature anisotropics in the CMB. 



4.5 QGP in the core of pulsars: 

Although I have focused the discussion of observable consequences of the QCD phase transition on the 
experimental signatures at SPS and RHIC, corresponding to a hot QGP, in section 3.7 I have presented 
arguments suggesting that the core of neutron stars may have a deconfined phase of quarks and gluons. 
In neutron stars the typical temperatures are of order 1 Mev while the baryochemical potential is of order 
~ Gev. Therefore a deconfined phase of quarks and gluons in this case correspons to a cold degenerate 

QGP m m n 



In a non-rotating neutron star, the boundary between the deconfined, mixed (hadron and QGP) and 
hadronic phases are fixed, but in a rotating neutron star (pulsar) these boundaries change as the rotational 
frequency of the star changes in time. In ref. |Q it was pointed out that since the compressibility of 
normal nuclear matter phase and the deconfined (almost free Fermi gas) of the QGP are different (the 
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Figure 10: HBT interferometry. 



incompressibility of normal nuclear matter is greater than that of an almost free Fermi gas of quarks) a 
structural change will occur upon a change in frequency. This has important consequences in the spin-up 



or spin-down stages of millisecond pulsars ||86| , 48 1. As a millisecond pulsar spins-down, its central density 



may rise above the critical density for the QGP phase transition in dense nuclear matter and the central 
core changes to a phase with a softer equation of state. In ref. [Q the argument is reversed to contemplate 
the spin-up stage of a X-ray pulsar in an low-mass X ray binary (LMXB's). In this case the X-ray pulsar 
accretes mass from its binary companion which is typically a low mass white dwarf, and spins-up. In this 
case if the density at the core falls below the critical density the, QGP at the core turns into the hadronic 
phase and the existing quark matter at the core is "spun-out". In this case, the frequency of these X-ray 
pulsars increases during accretion |Q. The authors of ref. |Q argue that for a range of frequencies the 
changes in the quark matter composition at the core will inhibit changes in the rotation frequency of the 
pulsar because of the increase in the moment of inertia. The net result of ref. is that these accreters will 
spend more time near these critical frequencies resulting in an anomalous distribution of frequencies at or 
near this frequency. A recent analysis of the oscillations of millisecond pulsars by the Rossi X-ray Timing 
Explorer (RXTE) |^ clearly shows that there is a frequency, ~ 300 Hz at which the pulsar distribution 
peaks. A detailed study of the time evolution of the moment of inertia and the rotational frequency for an 
LMXB was performed in ref. |Q. For the case in which the core of the pulsar has a deconfined QGP phase 
the authors show that the pulsar distribution has a spike at a frequency which is compatible with the RXTE 
data. The observed and theoretical distributions in frequency are shown in the figure ( pi] ) below. 

Glendenning and Weber Q argue that a degenerate quark matter core (cold and dense QGP) in the 
pulsars of some LMXB's of suitable masses can resist spin- up through the ongoing reduction of the quark 
matter cores in these accreting pulsars. As explained by these authors, a conversion from the quark to 
hadronic matter (the inverse situation as envisaged from a cooling QGP) manifests itself in an expansion of 
the star and a significant increase in the moment of inertia. The angular momentum added to the pulsar 
via accretion is used up by the star's expansion, inhibiting the spin-up until all the quark matter in the core 
has been transformed into the mixed (or hadronic) phase. 

These authors conclude that while there are several possibilities to explain the spike in the frequency 
distribution, this mechanism can contribute to explaining the anomalous frequency distribution. 

Other potential signals of a quark-hadron phase transition in neutron stars had been suggested, such as 
changes in the surface temperature as well as rotational mode instabilities . This wealth of potential 
observables can provide a definite astrophysical evidence of a deconfined phase of quarks and giuons in some 
of the most extreme environments in the present Universe. 

As mentioned above, there are fascinating novel color superconducting phases that are conjectured to arise 
in cold and dense QCD. The observational aspects of these novel phases of QCD are still being investigated 
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Figure 11: Frequency distribution of X-ray neutron stars. The spike in the calculated distribution is asso- 
ciated with the spinout of the quark matter phase and the corresponding growth of the moment of inertia. 
From reference p8[ 



and we refer the reader to the review articles in the literature for details. 

5 Back to the early Universe: summary 

The main goal of these lectures is to assess the observational possibilities of phase transitions in early Universe 
cosmology. I argued that the analysis of CMB anisotropics can reveal information on phase transitions that 
occurred during inflation. 

Phase transitions during the inflationary epoch modify the power spectrum of the primordial density 
fluctuations whose wavelengths cross the horizon during inflation and re-enter after recombination. Current 
theoretical models suggest that such phase transition occurred at a grand unified energy scale or perhaps 
within some supersymmetric theory but certainly beyond the current standard model and within a realm 
that neither theory nor experiment is on solid grounds. 

I have then narrowed the discussion down to the phase transitions predicted by the standard model of par- 
ticle physics because this model is on solid theoretical and experimental footing. However the observational 
consequences of phase transitions are argued to be rather indirect. 

The standard model predicts two phase transitions, one at the electroweak (EW) scale and the other at 
the QCD scale. The latter can actually be two phase transitions: a confinement-deconfinement PT between 
an almost free gas of quarks and gluons and hadrons and the other the chiral phase transition, both at a 
temperature T ~ 150 — 200 Mev. While originally it was conjectured that baryogenesis could be explained by 
a strong first order phase transition at the EW scale, current bounds from LEP for the Higgs mass seem to 



36 



rule out a strong first order phase transition. Furthermore, the magnitude of the CP violating parameters in 
the standard model, contained in the phase of the CKM matrix seem to be too small to lead to the observed 
baryon asymmetry. Thus, within the standard model the only phase transition that could lead to observable 
consequences are those of QCD. The QCD phase transition(s) occurred at a time t ^ 10^^ — 10^^ seconds 
after the Big Bang when the size of the horizon is dn 10 Km and the mass contained within the horizon 
is ~ 1 Mq. There are several possible consequences of the QCD phase transition that have been reviewed: 
i) Inhomogeneous nucleosynthesis, which could lead to a possible explanation of the new constraints for 
ilbh^ from the Boomerang/Maxima data analysis, solar mass primordial black holes, that could be part of 
the cold dark matter and strange quark nuggets or strangelets that can also be a component of CDM. All 
of these proposals are based on a first order phase transition and the details require a deep knowledge of 
the QCD parameters and equation of state. The lattice gauge theory program is providing reliable data 
on these issues but clearly there is a need for an experimental program to understand the feasibility and 
reliability of these observables. The experimental study of the EW phase transition itself is not feasible 
within our lifetime (certainly mine!) since energy densities ten orders of magnitude larger than that in 
nuclear matter are required. However the QCD phase transition(s) require energy densities a few times that 
of nuclear matter and are currently being studied by accelerators at CERN (SPS) and BNL (RHIC) with a 
forthcoming upgrade at CERN(LHC). 

Furthermore, after an excursion into the lives and deaths of stars I presented recent results that suggest 
that a deconfined phase of quarks can exist at the cores of neutron stars. 

I then summarized the program of relativistic heavy ion collisions that seeks to map the phase diagram 
of QCD and provided the recently reported analysis of "evidence" gathered at the SPS-CERN. While the 
interpretation of this evidence does not uniquely point to the discovery of a novel form of matter: the quark- 
gluon plasma (QGP), thus proving the confinement-deconfinement phase transition in QCD, taken together 
they represent a formidable body of "circumstantial evidence" in its favor. 

From the experimental data on J/\l/ suppression as well as photons and dileptons we will learn about 
the physics of the deconfined hot quark gluon plasma phase as well as the hot hadron phase. From data on 
strangeness production we can learn about the formation of strangelets during the QCD phase transition 
in the early Universe. From flow we can learn the Eos of QCD and if coherent pion domains are formed 
and measured we can then provide a more sound assessment on the possibility of primordial magnetic fields 
with net helicity seeded by the decay of these coherent pion domains. Thus the experimental program will 
undoubtedly lead to a firmer physical picture and a more solid basis for theoretical work. 

Furthermore, recent astrophysical observations of the frequency distribution of pulsars, which show an 
anomalous spike at ~ 300 Hz could in fact already be a telltale signal of a novel phase of cold and dense 
QCD in which deconfined quarks and gluons make up the core of pulsars. 

After SPS shuts down paving the way to the LHC, the torch is passed on to the Relativistic Heavy Ion 
Collider at BNL and the future LHC at CERN in which the ALICE program will study heavy ion reactions. 
Experiments at these colliders have the potential for understanding the details of the QCD phase transitions. 
Current and future astrophysical observations of X-ray spectra, timing and rotational properties of pulsars 
could confirm the possibility of deconfined quark matter at the core of these compact stars. Thus, indeed 
RHIC, LHC and pulsars are opening a window to the early Universe and giving a glimpse of its infancy when 
it was only 10~^ seconds young. 

Acknov^fledgements: The author thanks NSF for support through grants PHY-9605186, PHY-9988720 
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